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We calculate numerically the exact energy spectrum of the six dimensional problem of two in- 
teracting Bosons in a three-well optical lattice. The particles interact via a full Born-Oppenheimer 
potential which can be adapted to model the behavior of the s-wave scattering length at Feshbach 
resonances. By adjusting the parameters of the corresponding Bose-Hubbard (BH) Hamiltonian the 
deviation between the numerical energy spectrum and the BH spectrum is minimized. This defines 
the optimal BH parameter set which we compare to the standard parameters of the BH model. The 
range of validity of the BH model with these parameter sets is examined, and an improved analyti- 
cal prediction of the interaction parameter is introduced. Furthermore, an extended BH model and 
implications due to the energy dependence of the scattering length and couplings to higher Bloch 
bands at a Feshbach resonance are discussed. 



I. INTRODUCTION 

The Hubbard model has its origin in the description 
of electrons in solids [![. However, ultracold atoms in 
an optical lattice (OL) represent an almost perfect real- 
ization of this model with the additional advantage that 
many parameters such as the lattice depth and the in- 
teraction strength can be controlled and characteristics 
of the system can be observed with high accuracy 0. 
The Bose-Hubbard model (BHM) is a special form of the 
Hubbard model describing Bosonic particles in an OL. 
Despite its striking simplicity it is able to describe ex- 
citing phenomena such as the quantum phase transition 
between the Mott insulator and superfluid phase which 
is determined by the ratio U/J of the interaction and 
the hopping parameter [1, 0] . If the system is quenched 
from one quantum phase to another, local many-body re- 
laxation effects appear even without thermal relaxation. 
They are explained by the coupling between neighboring 
lattice sites in the BH Hamiltonian ||. Another predic- 
tion of the BHM is the existence of repulsively bound 
atom pairs in an optical lattice Q. In contrast to a 
real solid, the optical lattice does not allow dissipation 
to phonons. As there is no resonant unbound state in 
the lattice, the repelling atoms arc unable to decay. Ev- 
idence for this effect was found in an experiment with 
87 Rb atoms Q. 

Not only large optical lattices but also few-well systems 
such as double and triple wells are frequently modeled by 
a BH Hamiltonian [1,1^, [lfj. These systems have poten- 
tially a large range of applications. The triple- well system 
was proposed to serve as a transistor, where the popula- 
tion of the middle well controls the tunneling of particles 
from the left to the right well [ll|. Double wells, effec- 
tively generated by two interfering optical lattices form- 
ing a supcrlattice, are promising candidates to implement 
one- and two-qubit quantum gates [HI, EH, Q3 ■ 

To link experimental data to theoretical predictions of 
the BHM and to finally realize practical applications such 
as quantum gates, the precise knowledge of the parame- 
ters of the BH Hamiltonian is, however, important. We 



want to clarify in this work under which conditions the 
BHM is applicable and how the BH parameters may be 
improved. 

In order to examine and improve the validity of the 
BHM we calculate the energy spectrum of two interacting 
atoms in three wells of an OL by the help of a numerical 
routine introduced in Ref. [l5[. The reduction from the 
many-body problem to only two particles is senseful, be- 
cause the BHM is often used to describe dilute quantum 
gases with one or two atoms per lattice site. Also, the 
approximation of the OL by three wells does preserve all 
important microscopical features of the BHM of an OL. 
The BH Hamiltonian of the triple well allows for tunnel- 
ing of particles between the sites with amplitude J, for 
different onsite energies assigned to each of the sites 
i = 0, ±1 and of course for an onsite interaction with 
strength U. The spectrum of the BH Hamiltonian of the 
triple well consists of six eigenenergies which depend on 
the parameter set V = { J, U, eo, e±i}. Although the pa- 
rameters are predicted within the approximation of the 
BHM from the lattice depth, the confining potential, and 
the s-wave scattering length a sc , one can define an opti- 
mal set Vopt, which minimizes the deviation between the 
BH spectrum and the one obtained numerically. 

In this work we identify mainly two sources which can 
lead to large discrepancies between the predicted param- 
eters of the BHM and the optimal parameter set. First, 
effects of the additional potential which confines the par- 
ticles in three wells are clearly visible in shallow lattices. 
Albeit this is an artefact of the restriction to three wells 
it is, of course, important for the theoretical description 
of double and triple wells by the BHM. The second and 
more severe discrepancy emerges if the scattering length 
a sc is not significantly small compared to the character- 
istic trap length aho- As the Wannier basis used in the 
BHM does not reflect interaction it leads to a wrong pre- 
diction of the interaction parameter U for large scatter- 
ing lengths as they appear experimentally in the context 
of (magnetic) Feshbach resonances. In this regime the 
known analytical solution for two particles in a single 
harmonic well interacting via a point-like pscudopotcn- 
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tial is often used to describe effects of the interaction 
EE EM EH- If no perturbation theory is implied, 
the harmonic approximation is, however, limited to very 
deep lattices. In this work a simple-to-calculate correc- 
tion factor to the harmonic approximation is introduced 
that extends the validity to shallower lattices and im- 
proves the prediction even for deep lattices. 

Not only for strong interaction the standard BHM faces 
a potential weakness. For shallow lattices particles in 
neighboring wells are able to interact and hopping of 
particles to next-to-nearest neighbors can become impor- 
tant. One can, however, account for this by a so-called 
extended Bosc-Hubbard model. We examine in which 
regimes this extended model is necessary and useful. 

The paper is outlined as follows: In Sec. [IT] we describe 
the Hamiltonian of the triple well and we give a short 
review of the theoretical description of ultracold atoms 
in optical lattices described by the BHM. In Sec. [Hi] nu- 
merical methods and different approaches to determine 
and approximate the parameters of the BHM, especially 
the interaction parameter, are presented. Thereafter 
in Sec. IIVI the results of the full calculations and var- 
ious approximations are compared. Besides the values 
of the BH parameters and the resulting energy spectra, 
the influence of an energy-dependent scattering length 
is considered. Furthermore, the behavior at the res- 
onance of the scattering length including couplings to 
higher Bloch bands is discussed. Finally, we study ef- 
fects of next-neighbor interaction and hopping to non- 
neighboring wells. Conclusions are made in Sec. fVl 



II. SYSTEM 

A. Hamiltonian 

The two particles in the triple well are excmplarily 
represented by two bosonic Li atoms in the electronic 
a 3 S+ state. The Hamiltonian of this system in absolute 
coordinates is given by 

H [u n(ri,r 2 ) = -|- (v? + V|) + V^Qn - fi|) 
2m V / 

+ V tTap {n) + V tlap {f 2 ) (1) 

where Vt rap is the trapping potential and Vint de- 
scribes the 7 Li- 7 Li-interaction by the adiabatic Born- 
Oppenheimer (BO) potential of the a 3 Sj state. This 
electronic state relates to a sample of spin-polarized 
atoms. The corresponding BO potential has the advan- 
tage that it supports only few bound states which re- 
duces the number of basis functions needed to describe 
the wavefunction in the numerical calculation. 

For an ideal (isotropic) OL the trapping potential is 
given as 

Vb L (r)= E ^osin 2 (fc £) (2) 

£=z,y,z 



where Vq is the depth of the lattice and a = n/ko is 
its periodicity. The wavelength A = 2Tr/k Q is fixed to 
1 /xm in this work, which is in accordance to experiments 
with Li atoms in optical lattices [20| . The harmonic ap- 
proximation of the OL potential V(r) = ^moj 2 r 2 with 
u) 2 = 2Vofco/m defines the characteristic energy ftio and 
the characteristic trap length a\ lo ~ y / 2h/mcu which is a 
measure of the extent of the ground state solution in the 
harmonic approximation of the OL. Another character- 
istic unit of the OL is the recoil energy E r — h 2 k^/2m. 

To model a triple well, the OL potential is expanded to 
the twenty-second order in x direction and second order 
(i.e. harmonic approximation) in y and z directions. This 
leads to a new potential V22(r) consisting of three lattice 
sites which have in x-direction almost exactly the form of 
the OL potential (see Fig.[T]). The difference between the 
OL potential Vol and V22 can be regarded as an extra 
potential V con { = V22— Vol which confines the particles in 
three sites of an infinite OL. The triple- well Hamiltonian 
iJtrip for the two atoms N = 1, 2 is thus given as 

£ (§; + v& + v<l?l+v^. 0) 

2V=1,2 ^ ' 

In a real experimental setup the confinement to three 
wells could be either due to a superlatticc or could be 
generated by the harmonic potential of an optical dipolc 
trap which is, however, usually much shallower than in 
our case where V con f(r) rises rapidly for \kox\ > 
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FIG. 1: (Color online) Visualization of the "hopping" -parameter 
J, the onsite energies e; (i = 0, ±1) and the interaction strength 
U. The potentials V22{r) (blue solid), Vol^) (black dashed) and 
Vconf (brown dash-dotted) are depicted for y = z = 0. 



B. Bose-Hubbard model 

The periodicity of the OL potential gives rise to band 
gaps in the dispersion relation E(k). In contrast to a 
real solid, both the dispersion relation and the Bloch so- 

lutions of the Hamiltonian -ffoL = ~^ Vol are known 
analytically. The extents of the first four Bloch bands of 
a one-dimensional OL are depicted as a function of the 
lattice depth Vq in Fig. O For increasing lattice depths 
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the band widths shrink and the gaps between the bands 
increase. Especially, the gap between the first and the 
second band prevents ultracold Bosons from occupying 
others than the lowest Bloch band. The energy of states 
in this band can be roughly approximated by the ground 
state energy ftw/2 of the harmonic approximation. How- 
ever, even for Vo ~S> fico there is a constant energy offset 
of E r /A between the harmonic approximation and the 
exact ground state energy. This offset is a result of the 
anharmonicity of the lattice sites. It can be explained 
by a perturbation of the harmonic ground-state energy 
by the next-to-leading order expansion (~ x 4 ) of the OL 
potential 21 1. 



V [£,] 




v„ PM 

FIG. 2: (Color online) Extents of the lowest four energy bands 
of the one-dimensional optical lattice as a function of the lattice 
depth Vq. The energy of the first band can be approximated by 
the harmonic ground-state energy 5^ (dashed line). A better 
approximation includes a constant energy offset of E r /4 (see the 
text). 

The BHM makes use of the basis of Wannier functions, 
constructed with the aid of the known Bloch solutions of 
the lowest band to approximate the triple-well Hamilto- 
nian H tr i p . In the case of an infinite OL with no interac- 
tion this orthogonal basis spans completely the Hilbcrt 
space of the first Bloch band. In contrast to the Bloch so- 
lutions, each Wannier function Wj(r) is localized at one 
lattice site j. Wannier functions can be constructed in 
many different ways depending on the choice of phases 
of the Bloch solutions. In accordance to the widely used 
convention of the BHM we consider the form specified 
by Kohn in [22| which leads to maximally localized, real 
Wannier functions. 

To describe systems with a finite number of wells, the 
amount of Bloch solutions is usually restricted to those 
which fulfill periodic boundary conditions. For the triple- 
well potential used in this work, the boundaries may be 
set to koXB = ^f- or, in order to probe the outer potential 
walls, to koXB = 27r. These conditions change slightly the 
shape of the Wannier functions. 

We adapt the basis to the triple- well potential by using 



in y- and z-direction the appropriate harmonic ground- 
state solution ho instead of the Wannier function for a 
one-well lattice. The adapted Wannier functions read 



Wj(f) = Wj (x)h (y)h (z). 



(4) 



A detailed description of Wannier functions can be found 
in Refs. [Iim. 

Let bj (bj ) be the bosonic operator of creation (annihi- 
lation) of a particle with the Wannier function Wj(r) = 
(f\wj), then Atrip (Eq.|3J) is written in the Wannier basis 
as 



Wan 



i W *\ f- + ^OL + Vconf \Wj) b\bj+ 
^ (w k \(wi\ V int \w m )\w 



b\b\b m b n , 



(5) 



with all indices running over —1,0, +1. 

This is already the first approximation of the BHM as 
the Wannier basis is only complete for the first Bloch 
band if V int = and V con [ = 0. The BHM further simpli- 
fies .ffwan by making the following assumptions: (i) The 
overlap between Wannier functions on different lattice 
sites is small. This reduces the i, j-summation to diago- 
nal elements i = j and neighboring lattice sites (ii) 
The confinement potential V con f varies slowly and does 
not couple Wannier functions at different lattice sites, 
i. e. (wi \ V con f \wj) = for i ^ j. (iii) The interaction 
is short-ranged and only takes place between particles in 
the same well. 

Applying these assumptions, the final BH Hamiltonian 
takes the form 

Hbh = -J + f E W b *i + E e ^k (6) 

<i,j> i i 



with the hopping parameter J = — 
Vol the onsite energies e 4 = {Wi\ + Vol + 

^conf \wi) (i = 0, ±1) and the interaction parameter U = 
(wq \ (wq \ Mnt I wo) |u>o). Even though the BH Hamilto- 
nian has a simple form, there does not exist a general 
solution for large lattices with many particles. However, 
for two particles in a triple well the BH Hamiltonian is 
easily diagonalizable. 

The integrations occuring in the definition of the BH 
parameters are performed consistently within the peri- 
odic boundary conditions of x. The Bloch functions 
which are used to construct the Wannier functions are 
also normalized to one within the boundaries. Since the 
trapping potential is symmetric, e_i = e+ 1 holds. For a 
visualization of the BH parameters see Fig. [TJ 

A further important approximation, which is very com- 
mon in the physics of ultracold atomic gases, is the re- 
placement of the full BO interaction potential by a point- 
like pseudopotential [24[ 

T . / —*\ 47rft 2 a sc ./,-, ._. d 
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that depends only on the s-wave scattering length a sc . 
This simplifies the integration for the interaction param- 
eter drastically to 

[f=^/dV W f)f, (8) 
m J 

The formulations of J,U,£q, and e±i described above 
are in the following denoted as the "BH representation" 
of the parameters and will be indicated with the super- 
script "BH". 

It should be noted that not only the Wannier ba- 
sis leads to the form of the BH Hamiltonian (Eq. [6]). 
Any other nearly orthogonal basis of functions local- 
ized in one of the wells leads to the same form, 
with only different values of the parameters. The 
usual basis of the BH Hamiltonian of two parti- 
cles in a triple well is the occupation-number ba- 
sis {|002) , |011) , 1 101) , |020) , |110) , |200)} which speci- 
fies the amount of particles in the left, middle and right 
well. In this respect it is even possible to think of a 
two-particle instead of a one-particle basis to arrive at 
the BHM. As a result, the BH Hamiltonian is principally 
able to describe strongly interacting particles without the 
need to expand the full Hamiltonian (TT|) in a large basis 
of one-particle Wannier functions of several Bloch bands. 
It is hence natural to imagine that there is a basis for the 
BHM which is far better adapted to an OL with inter- 
acting particles in the presence of a confining potential 
V^ on f than the Wannier basis. Although this optimal ba- 
sis is unknown, it allows us to vary the BH parameters in 
order to adapt the model to the full numerical solution 
as is described below. 

III. METHOD 

A. Numerical calculation in the triple well 

To perform the full numerical calculation of the 
cigenencrgies and eigenfunctions of Hamiltonian JTtrip 
the approach described in [l5| is adopted. Briefly, in this 
method the lattice potential is expanded into a Taylor 
series and the resulting Hamiltonian is written as a sum 
of center-of-mass (COM) Hamiltonian i/coM, relative- 
motion (REL) Hamiltonian Hrel and coupling terms W. 
The eigenfunctions and eigenvalues of Hqom and i^REL 
are calculated separately in a basis of spherical harmonics 
Y^iO, (f>) times B splines B a (r). The position of knots of 
the B splines in REL can be adjusted to adequately de- 
scribe both the long range behavior in the trap (~ 10 a>o) 
and the short range behavior due to the molecular (inter- 
atomic) interaction (~ 65ao). The full solution including 
the coupling W is calculated from the eigenfunctions of 
ffcOM and Hrel by the configuration-interaction (CI) 
method which is also known as exact diagonalization. 
The code is adapted to the symmetry group Z?2h of a 
cubic three-dimensional OL, which simplifies the calcu- 
lation and allows to reduce the REL basis to functions 



which arc symmetric under inversion and describe there- 
fore Boso nic p articles. 

In Ref. [la ], a Taylor expansion up to the sixth order 
of the OL is used to examine the two-particle interac- 
tion in a single well. The ground state was described 
by a basis which consisted of angular components up to 
1 = 3. The triple well considered in the present work 
demands much higher angular components owing to its 
high degree of anisotropy. We need an expansion up to 
I = 32 and about 70 B splines in COM and an expan- 
sion up to I = 30 and 130 B splines in REL to obtain 
converged solutions for lattice depths up to Vq = 12E r . 
Altogether, about 16 000 basis functions are used to de- 
scribe the REL eigensolutions and 1 1 000 to describe the 
COM eigensolutions. From these solutions about 10 000 
configurations form the basis of the CI calculations. 

To tune the s-wave scattering length a sc to values in 
the range a sc € (— oo, +oo) the inner wall of the BO 
potential is slightly shifted as suggested in Ref. [15[ . This 
procedure can move the least bound state close to the 
dissociation threshold or leads to the creation of a new 
bound state. Simultaneously the scattering length grows 
to +oo or -co, respectively, which offers a knob to vary 
a sc in a controlled way. 

B. Numerical Determination of the optimal BH 
parameters 

The numerical determination of the optimal BH pa- 
rameter set is done by fitting the parameters V = 
( J, eo, e±i, U) such that an optimal agreement between 
the eigenenergies of the BH Hamiltonian and the full nu- 
merical solution is achieved. 

In more detail, let 

#trlp|**}=-Ei|*<) (9) 

be the stationary solution of the triple-well Hamiltonian 
and 

h£k\$Y) = e? (io) 

be the stationary solution of the BH Hamiltonian Hg H = 
-Hbh with parameter set V . Due to the restriction to one 
basis function per particle and lattice site, the second 
eigenvalue equation has only six solutions with energies 

cV cV 

c 1 , . . . , c 6 . 

With the aid of the eigenenergies, we obtain the opti- 
mal set of BH parameters "P opt by minimizing the mea- 
sure 

i— 1 x 7 

However, a fit of all four parameters at once is problem- 
atic. Since only six eigenenergies exist, a local minimum 
of f(V) may be found that leads to a non-optimal pa- 
rameter set V. On the other hand, it is inherent to the 
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BHM, that all effects due to interaction arc contained 
in the interaction parameter U leaving the other param- 
eters unaffected. Therefore, we determine the optimal 
parameter set in a two-step process: (i) The parameters 
J> £ 0j e ±i a re fitted to match the energies of a numerical 
calculation without interaction, (ii) The parameter U is 
fitted to match the energies of a numerical calculation 
with interaction, keeping the other BH parameters fixed 
to the values of the previous fit. 



C. Estimation of the interaction parameter U 

As already mentioned, the Wannier functions do not 
reflect effects of interaction. Being a one-particle basis, 
the Wannier basis is only complete for non-interacting 
particles. Thus, one cannot assume that in strong inter- 
acting regimes the BH representation of the interaction 
parameter given in Eq. ([5]) is applicable and one should 
estimate the energy offset caused by interaction differ- 
ently. 

Fortunately, an analytical solution is known for a 
system of two identical particles in an isotropic or 
anisotropic harmonic t rap interacting only via the pseu- 
dopotential of Eq. ([7]) 0, [2{| . While for no interaction 
(a sc = 0) the eigenenergies in the isotropic case are given 
as E n _i = 7iw(§ + 2n + I), a non- vanishing scattering 
length leads to a shift of all eigenenergies with 1 = 0. 
The dependence of the energies -E ra ,z=o on the scattering 
length is shown in Fig. [3] and is given by the solution of 
the implicit equation 



1 



tan 



Oho 



/7re 7r\ 

(y + i) 



n £(j+j) 

4^(1 + 1) 



(12) 



where e = E/hu) and T denotes the Gamma function 
[26| . Additionally, for positive scattering lengths one 
bound eigensolution appears (lowest branch). For ultra- 
cold atoms only the first two states with energies lower 
than £?o,i = f should be occupied. This corresponds 
to the reduction to the first Bloch band in the OL. 

The known energy spectrum offers the possibility to es- 
timate the energy offset due to interaction as a function 
of the scattering length. As shown in Fig. [3l for posi- 
tive scattering lengths one has the choice between two 
branches; one with a negative offset AE® m to the lower 
branch and one with a positive offset AE^ arm to the up- 
per branch. For negative a sc only the latter branch exists 
and Ai?^ arm is negative. 

Since the interaction of the BHM has purely onsite 
character, one can expect that these energy offsets pro- 
vide a good first approximation [/ harm — AE®'^ Tni of the 
interaction parameter for strong interaction. Depending 
on which offset is chosen lower or higher parts of the en- 
ergy spectrum may be modeled. The confinement in the 
harmonic trap is, however, stronger than in a sinusoidal 
OL. Thus, one can expect that the harmonic approxi- 
mation tends to overestimate the strength of the onsite 
interaction in the OL. 
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FIG. 3: (Color online) Analytical energy spectrum of relative- 
motion £ = states of two particles interacting by the potential V pB 
(Eq. in a harmonic trap as a function of the scattering length. 
The spectrum defines the energy offsets AE? and AE} rela- 

1 °^ harm harm 

tive to the ground-state energy in the absence of interaction, ^Hu>. 



Despite the weakness of the Wannier basis for strong 
interaction it can offer crucial information on anharmonic 
effects. Its completeness for Vint = and y con f = al- 
lows an almost exact prediction of the behavior of U for 
a sc — > 0. (To first order, V con { has no influence on U.) 
This observation can be used to define a correction fac- 
tor A to the harmonic approximation which accounts for 
anharmonic effects. The factor is defined by demanding 
that A ■ AE^ mm has the same linear behavior as £/ BH for 
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(13) 



Remarkably, the factor A, which is shown in Fig. U for 
different lattice depths Vo, approaches even for very deep 
lattices only slowly unity. Thus it is needed to correct 
the harmonic approximation even for very deep lattices. 
For shallow lattices the factor rises above one and goes 
for Vq — ► to +oo. This is due to the fact that periodic 
boundary conditions are implied. Thus, the normaliza- 
tion of the Wannier functions within the boundaries pre- 
vents the integral in the BH representation of U to vanish 
for Vo — > while a^o/huo goes to infinity. 

The validity regime of the corrected harmonic approx- 
imation [/ corr = .4[/ harm depends on the applicability 
of the pseudopotential approximation with an energy- 
independent scattering length used in [/ harm and the va- 
lidity range of A. The latter can be roughly predicted by 
considering that the correction factor accounts for the 
different confinement in the harmonic trap and the real 
OL for a certain ground state energy E of noninteract- 
ing particles. In order for A to be valid for a sc ^ 0, the 
deviation of the energy (ss U) should be small compared 
to the depth of the lattice Vo, which simply means 



\U\ w \U co "\ < V . 



(14) 
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stron g co nfinement and strong interaction in harmonic 
traps [H HIOH . 

To quantify the impact of the energy dependence in- 
dependently of anharmonic effects we consider two 7 Li 
atoms in a harmonic trap and determine numerically the 
ground-state energy. Plugging this energy into the right- 
hand side of Eq. (p~2|) defines the "optimal" scattering 
length a opt for the description of the interaction by the 
pscudopotcntial. In its range of validity the effective scat- 
tering length is equal to the optimal scattering length. 

The numerical calculation with the full BO interaction 
in an OL makes it possible to compare the effect of an 
energy-dependent scattering length to effects caused by 
the anharmonicity of the OL and the incompleteness of 
the Wannier basis. 



FIG. 4: (Color online) The correction factor A as a function of the 
lattice depth V . Inset: Comparison of U Bn , (7 harm , and U COIT 
for Vo = fiui as a function of the scattering length. In accordance 
with the definition of A the behavior of f/ BH and £/ corr matches 
for small a sc . 



This is to some extent counter-intuitive as for a fixed 
scattering length a sc an increase of Vq is accompanied by 
an increase of |a sc /aho|- This leads to a breakdown of the 
BH representation and to errors caused by the use of an 
energy-independent scattering length in [/ harm . Since at 
the same time the error of the correction factor A shrinks 
it is important to quantify and compare the different er- 
rors in order to understand when JJ co " is applicable. 



D. Sextic approximation 

An alternative way to look at effects of anharmonic- 
ity on the interaction is to consider an optical lattice not 
expanded to second or twenty-second but to the sixth or- 
der in x-direction. This sextic trap consists of one well, 
but includes anharmonic effects and coupling between 
relative and center-of-mass motion (l5j . The energy dif- 
ference of two atoms in a sextic well with and without 
interaction defines the interaction parameter [/ sext . 



E. Energy-dependent scattering length 



IV. RESULTS 

This section is devoted to a comparison of the ana- 
lytical approaches to describe interacting particles in a 
triple well with the full numerical results which deter- 
mine the optimal BH parameters. First, the interaction- 
independent parameters J, eo, e±i arc considered for a 
varying lattice depth Vq. Then the discussion is extended 
to interacting particles for both different lattice depths 
and different scattering lengths. 

Here a short overview of the abbreviations used in the 
following. 

BH Parameters in BH representation and 
eigenenergies of the BHM with these 
parameters. 

harm The same as "BH" but with [/ harm deter- 
mined in the harmonic trap, 
corr The same as "BH" but with the corrected 

rycorr t r^harm 

opt The optimal BH parameters and eigenen- 
ergies of the BHM with these parameters. 
This abbreviation is also used for the op- 
timal scattering length. 

sext Interaction parameter JJ scxt obtained 
from the sextic approximation of the OL. 

num The eigenenergies of the full numerical CI 
calculation. 



The scattering length used in the pseudopotential 
(Eq. [7]) is defined by the s wave of a trap-free scatter- 
ing process in the limit of zero collision energy The so- 
lution in the trap has, however, a different asymptotic 
behavior and a finite ground-state energy. Much work 
has been devoted to the question when a pseudopotential 
can be used in a trap and which scattering length should 
be chosen for energies greater than zero [H, |2g, H3, HU ■ 
It has been shown that the use of an energy-dependent 
effective scattering length can largely extend the range 
of validity of the pseudopotential approximation towards 



A. Onsite energies 

The optimal onsite energies for the middle well and the 
outer wells are shown in Fig. [5] On the first glance, they 
seem to be very well matched by the corresponding BH 
representations. However, by subtracting the value of the 
onsite energy for no confinement (wq\ l^o) m t ne 

inset of Fig. [5] the impact of V con { on the BH representa- 
tion and the optimal values of the onsite energies can be 
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examined. Generally, the shallower the lattices the more 
the optimal onsite energies of both wells are lifted by the 
confinement. Only for Vb > hui this trend is reflected by 
the BH representation, while for shallower lattices bound- 
ary conditions at k$XB = 3tt/2 underestimate the impact 
of the confinement. On the other hand, for k$XB = 2ir 
(not shown in the graph) the boundary reaches deep in- 
side the walls of the confinement and the BH represen- 
tation extremely overestimates e±± for Vq < 1.7 hui and 
even eo for Vq < 0.8 hu> . 

A reason for the inaccuracy of the BH representation 
of eo and e±i is that only for a slowly varying confining 
potential V COIl {(f) the Wannier functions form a good 
basis. As V con f(f ) varies strongly for \kox\ > ^ the 
lattice must be considerably deep so that the particles 
do not "see" the form of the boundary and the distance 
Xb of the boundary conditions must be sufficiently small 
not to reach inside the region of the strong increase of the 
confining potential (see Fig. [l}. For boundary conditions 
at k xs = 3ir/2 and lattice depths Vq > 1.5 hu> both 
the values of the onsite energies and their two orders of 
magnitude smaller differences are well described by the 
BH representation. It can be assumed that for larger 
lattices with smooth harmonic trapping potentials, the 
BHM converges faster. 
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FIG. 5: (Color online) Onsite energies as a function of the lattice 
depth Vb : the optimal values together with the BH representation 
for the middle well (eo) and the outer wells (e±i) f° r koxs = 3tt/2. 
Inset: In order to examine the impact of the confining potential 

-2 

V con f the value of 8ei = ti — (wq\ J— + Vol \il>o) is depicted, where 
the value of the onsite energy for no confinement is subtracted. 
Clearly, for Vq > 1.5 hui both results are converged. 



B. Hopping parameter 

As in the case of the onsite energies, the confining 
potential can cause also problems for the prediction of 
the hopping parameter. Fig. [5] (a) shows the optimal 
hopping parameter and its BH representation for differ- 



ent boundary conditions k^xs — 37r/2,27r and oo. One 
can see that the choice of the boundary conditions in- 
fluences considerably the BH representation for shallow 
lattices. It seems that in this special case boundaries 
at 2tt or infinity produce the best prediction. However, 
similarly to the case of the onsite energies all BH rep- 
resentations predict the correct hopping parameter for 
Vq > 1.5 hui. The reason is analogous. The BHM as- 
sumes that V'conf is approximately constant within one 
well. Therefore, due to the orthogonality of the Wan- 
nier functions, (wq\ V con { \wi) = holds. This is essen- 
tial for the definition of the hopping parameter, where 
the influence of the confinement potential is neglected. 
However, as described in the last paragraph, for shallow 
lattices the form of the boundary becomes important and 
strongly varying parts of the potential are probed. One 
can expect an error of the BH representation for J of 
the size A J = (wq\ V corl { \wi) which is negligible only for 
sufficiently deep lattices. 



Energy spectrum for no interaction 

Since the BH representation of both the hopping 
and the onsite-energy parameter is converged for Vq > 
1.5 hui, the BHM predicts the correct eigenenergies in this 
regime. The mean deviation from the numerical energies 
for Vo = 1.7 Hu is 0.25£V or 0.01% respectively Fig. ® 
(b) shows the spectrum of the full numerical calculation. 
In order to compensate for an increase of the energies 
of the order of 2eo and the shrinking width of the first 
Bloch band of the order of 4 J as Vb is increased, rescaled 
energies r\i with 



rn(y ) 



E, 



2f BH 



4J BH 



(15) 



are shown. For Vq > 1.5 hu> the spectrum of the BH 
Hamiltonian in the BH representation is in very good 
agreement with the spectrum obtained from a full nu- 
merical calculation. With the optimal parameter set the 
BHM is able to predict the numerical eigenenergies al- 
ready for Vq > hco with high accuracy (mean relative 
error in energies w 0.001%). For shallower lattices nu- 
merical and optimal eigenenergies do not match, which 
means that the BHM is principally unable to reproduce 
the eigenenergies of the triple well for Vq < huj. Two 
differences are most obvious: (i) Both the BHM in the 
BH representation and with the optimal BH parameters 
underestimate significantly the value of the ground-state 
energy while especially the BH representation tends to 
overestimate all other eigenenergies in the limit of van- 
ishing lattice depth, (ii) For (7 = the BHM predicts an 
energy difference of e\ — eo between the third and fourth 
energy level. In the BH representation this difference 
goes to zero for Vq — ► (see Fig. [5]) and the levels be- 
come degenerate. Although this degeneracy is lifted by 
using the optimal parameter set, the energy difference is 
still largely underestimated. 
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0.5 1.0 1.5 

FIG. 6: (Color online) (a) Hopping parameter J as a function 
of the lattice depth Vb : the optimal value together with the BH 
representation for koXB = §T, 2tt, oo. (b) The six eigenenergies in 
the rescaled form (see Eq. I|15|l ) for a BC = and U = as a function 
of the lattice depth Vb- The BH energies are obtained for onsite 
energies with boundaries at koxs = 3tt/2 and a hopping parameter 
with boundaries at kgxs = 2-7T. 

These disagreements can be caused by the dominance 
of the confinement potential over the OL potential for 
small lattice depths. Also hopping to non-neighboring 
wells can become important as shown in Sec. IIVDI As 
already mentioned, one can assume that the BHM con- 
verges better for larger lattices with a shallow harmonic 
external confinement. On the other hand, the conver- 
gence behavior for, e.g., double wells as they appear in 
superlattices should be similar or even worse compared 
to the case of the triple well. 



C. Interaction parameter 

The interaction parameter U depends on both the scat- 
tering length a sc and the lattice depth Vq that determines 
how strong two particles are confined in one well. First, 
the dependence on the lattice depth for the case of a rel- 



0.5 1.0 1.5 

V [fuo] 

FIG. 7: (Color online) (a) The Interaction parameter as a function 
of the lattice depth Vb in the weak interacting regime with a BC ~ 
— 180ao. This results in 0.02 < \a BC /a^ \ < 0.08 in this graph, 
(b) The six eigenenergies in the rescaled form (see Eq. ||16| |) as a 
function of the lattice depth Vq. 



atively small and a relatively large negative scattering 
length will be discussed. 

Fig. [7] (a) shows a plot of the optimal interaction pa- 
rameter U opt together with several predictions of the in- 
teraction energy by U Bii , U halm and U co " . The scatter- 
ing length is fixed to about — 180ao- For this value, the 
ratio |a sc /aho| which increases with the lattice depth Vb 
reaches up to only 0.08 in this graph. 

Within the BHM the interaction parameter J7 BH is in- 
dependent of the confining potential V^onf- Accordingly, 
the main source of error is not that the Wannicr functions 
are not adapted to the confining potential, but that they 
do not reflect effects of interaction. Thus, for weak in- 
teraction the BH representation predicts U opt quite well 
already for shallow lattices with Vq ~ hco. However, an 
increasing value of Vq amplifies the confinement and thus 
the ratio |a S c/&ho|- Remarkably, this causes a small error 
of C/ BH which is already visible for |a sc /aho| < 0.08 and 
reduces the validity regime of t/ BH to scattering lengths 
which are at least one order of magnitude smaller than 
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FIG. 8: (Color online) (a) The Interaction parameter as a lunction 
ol the lattice depth Vo in the strong interacting regime with a EC ~ 
— 4600ao. This results in 0.5 < |a sc /ai 10 | < 2.2 in this graph, 
(b) The six eigenenergies in the rescaled form (see Eq. Jl6j) as a 
function of the lattice depth Vq- 



the trap length dh - 

For weak interaction the energy-offset in the harmonic 
well [/ harm can predict U opt only poorly compared to 
[7 BH . On the other hand, the corrected analytical ap- 
proach JJ C °" introduced in Sec. IIII CI is in excellent 
agreement with the optimal interaction parameter for 
Vq > hoj. This is also reflected by the energy spectrum 
shown in Fig. [7] (b). 

In Fig. [JJ (b) the upper three states describe particles 
in different wells with energies almost independent of the 
interaction. The first three states, on the other hand, de- 
scribe attractively interacting particles in the same well. 
Within the BHM their energies are shifted in the order of 
U. Thus, the rescaling of the energies (Eq.[T3]) is changed 
to 



BH representation and the harmonic approximation un- 
derestimate the energies more and more as the lattice 
depth and with this |a sc /aho| increases. For lattice depths 
smaller than huj the errors from the non-interacting case 
are inherited, e. g. the ground-state energy and the split- 
ting between the third and the fourth level is underesti- 
mated by all models. 

In Fig. [5] (a) the interaction parameter is plotted for a 
large negative scattering length of about — 4600a which 
means, that the absolute value of the scattering length is 
up to two times larger than the trap length a\ w . As the 
magnitude of a sc is of the order of aho, the BH represen- 
tation of U differs drastically from the real interaction 
parameter. In this regime the linear dependence of J7 BH 
on the scattering length leads to significantly wrong re- 
sults (see the inset in Fig. H}. In this case, the energy 
offset fj harm provides a far better approximation, but it 
results in an almost constant underestimation of the opti- 
mal interaction strength as is evident from Fig. [5] (a). The 
corrected interaction parameter, on the other hand, gets 
closer to U° pt for deep lattices. The graph supports that 
Eq. (fT4")l indeed reflects the validity regime of A. As U is 
approximately proportional to Vb/fiw the ratio Vo/|J7 corr | 
shrinks with an increasing lattice depth like {hw)^ 1 and 
U com ' approaches the optimal interaction parameter. The 
same is again reflected by the energy spectrum shown in 
Fig. 0(b). Only U c °" is able to reproduce the numerical 
eigenenergies of the three lowest states in deep lattices, 
while ;/ harm generates an almost constant underestima- 
tion and ?7 BH a rapidly increasing underestimation of the 
numerical eigenenergies of these states. 



1. Sextic approximation 

Figs. [Jj (a) and [8] (a) show additionally the estimate 
of the interaction parameter using JJ sext , the interaction 
energy obtained in the sextic trap. It agrees precisely 
with the optimal parameter U opt independently of the 
strength of interaction provided that Vo > 1.5 hoj. Also 
for shallower lattices, it is close to the optimal value. 
Thus, the examination of the interaction in a sextic trap 
[l5| can already cover all anharmonic features of the on- 
site interaction for sufficiently deep lattices. This implies 
that one can treat the anharmonic form of a lattice site as 
a perturbation to the harmonic approximation to study 
interactions in OLs more realistically. For example, in 
[l9l ] this has been done by treating the difference between 
the OL potential and the harmonic approximation as a 
perturbing potential. The effects of the anharmonicity 
within a single site of an optical lattice have also been 
investigated non-perturbatively in [30| . 



VIM) 



E, 



2f BH 



4J BH + \U C 



(16) 



2. Energy-dependent scattering length 



The energies of the lowest three states are well pre- 
dicted by using JJ colr for Vq > 1.2 fku while both the 



One could expect that the use of the optimal energy- 
dependent scattering length a op t instead of a sc can im- 
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prove the value of [7 BH (see Eq. (|8])) significantly. This 
is, however, not the case, as can be clearly seen in Fig. [7] 
(a) and Fig. [8] (a). In the validity regime of the BH rep- 
resentation (|a sc | <C aho) the energy dependence of the 
scattering length seems not to be of importance. Also 
for strong interaction the offset between C/ BH when used 
with a sc or a opt is small compared to the offset between 
any of the Z7 BH and U° pt . In other words, the effect of 
the energy dependence of the scattering length in the OL 
is generally small compared to the error due to the use 
of the interaction-independent Wannier basis. 

To quantify the effect of the energy-dependence, we 
use 



AU 



\a sc ) - t/ harm (a opt ) 



C/ harm («opt) 



(17) 



AU. Therefore, the total error of U COT1 ' is generally much 
smaller than the error of £/ BH or [/ harm . 



TABLE I: The error AU of the use of an energy indepen- 
dent scattering length together with errors of U BH , [/ harm and 
rycorr re i a tj ve to U° p for various lattice depths and scattering 
lengths. 



V /E r 


11.5 


11.5 


25.2 


64.7 


Osc/aho 


-0.08 


-2.01 


-2.44 


-3.09 


U COII /V 


-0.05 


-0.42 


-0.30 


-0.20 


AU 


-0.58% 


-1.5% 


-1.8% 


-2.2% 


Error of t/ BH 


-3.18% 


-190% 


-243% 


-320% 


Error of C/ lmrm 


-8.53% 


-7.0% 


-5.1% 


-4.2% 


Error of C7 corr 


0.06% 


1.5% 


0.04% 


-1.1% 



the error of the interaction energy in the harmonic well, 
if a sc is used instead of a opt . This estimates the error of 
neglecting the energy dependence of the scattering length 
independently of anharmonic effects. For a sufficiently 
large value of |a sc /aho| the deviation of [/ harm from the 
linear behavior of L7 BH (see inset of Fig. 0]) starts to be 
significant. Therefore, the size of a sc relative to a op t does 
not directly show the impact of the energy dependence. 

In Tab. U besides AU the errors of [/ BH , U haim and 
jycon- com p arcc | to the optimal value are listed for vari- 
ous scattering lengths and lattice depths. For very deep 
lattices with Vq > 12E r where a full numerical calcula- 
tion in the triple well is too laborious we use the sextic 
approximation U scxt instead of the optimal interaction 
parameter, as it exactly reproduces U opt for Vo > 10E r . 

The value of AU in Tab. U shows to depend mainly 
on the ratio a sc /aho- However, even for a large ratio 
it amounts up to only a few percent. Also the error of 
f/ BH depends on that ratio, but grows much faster with 
|«sc/flho|- This reflects again the complete breakdown 
of the BH representation for scattering lengths with a 
magnitude of the order of the trap length. 

The error of t/ harm decreases, though quite slowly, with 
growing lattice depth Vq. It comprises the error of an 
energy-independent scattering length (AU) and of ne- 
glecting anharmonic effects. The latter effect can be es- 
timated by subtracting the value of AU from the error 
of [/ harm . Only for the largest depth considered in Tab.[fl 
(Vq = 6A.7E r ) effects of AU are bigger than anharmonic 
effects. For a fixed lattice depth Vo = 11.5E r and a de- 
crease of the scattering length from — 0.08ah o to — 2.01ah o 
anharmonic effects seem to become less important. This 
can be explained by the stronger confinement in the wells 
due to the attractive interaction. 

The error of U corl ' results from the error of the correc- 
tion factor A and from the error AU. The error of A 
can again be estimated by subtracting AU from the er- 
ror of U colT . This shows that the error of A depends, as 
expected, on the ratio U corl /Vq. In all considered cases 
it is found to be smaller than the error of £/ harm caused 
by the anharmonicity and possesses the opposite sign of 



3. Interaction energy at a resonance of the scattering length 

Until now the two cases of very small and very large 
interaction were discussed. To examine in more detail the 
regime of validity for [/ BH , [/ harm an d [/con- we congjjgj- 
in Fig. [5] the behavior of the interaction parameter at 
a resonance (dho/flsc — > 0) for a lattice depth of Vq — 
X.lhio. At this depth the BH representation of the other 
parameters (J, eo and e±i) is close to their optimal values. 

Since at a resonance the energy regime of higher Bloch 
bands can be reached the notion of Bloch bands should 
be reconsidered in the presence of the new interaction 
parameters [/ harm and U c °" . For strong interaction one 
normally speaks of a coupling to higher Bloch bands, 
which means, that the true eigenfunctions can only be re- 
produced by an expansion into Wannier functions of sev- 
eral Bloch bands. The new interaction parameters are, 
however, not formulated as matrix elements in the Wan- 
nier basis. The resulting BH Hamiltonian is therefore 
intrinsiccly a multi-band Hamiltonian. As shown in the 
last paragraphs the use of an optimal parameter set can 
reproduce the numerical energy spectrum very well. One 
can therefore describe multi-band physics surprisingly 
well within the simple BHM. As we will show now this 
has, however, limitations for repelling particles where the 
BH states become degenerate with excited "molecular" 
states and noninteracting ones in higher Bloch bands. 

First of all, by varying the scattering length over the 
resonance (aho/«sc = 0) it becomes obvious, that the 
BH representation [/ BH is not valid in the proximity of 
the resonance where it diverges to ±oo (see Fig. |9] (a)). 
On the other hand, U c °" approximates well the optimal 
value U° pt for a large range of repelling (upper branch 
in Fig. [9] (a)) and attracting (lower branch in Fig. [9] (a)) 
interactions. 

In Fig. [9] (b) the lower branch is shown on an enlarged 
scale around the resonance. It shows that U COTT is valid 
until |a sc | ~ aho- For growing \U\ the correction factor 
A looses slowly its validity. Again, one can observe that 
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FIG. 9: (Color online) The Interaction parameter for different 
scattering lengths at ^ = 1.7. To observe the behavior at the 
resonance, the dependence on 2ta i s shown. 



the behavior of U° pt is exactly followed by the sextic ap- 
proximation while the harmonic approximation generally 
overestimates the magnitude of U opt . 

On the upper branch in Fig. [9] near the resonance the 
energy of bound particles in higher Bloch bands can cross 
the BH states. To understand how this happens, we have 
sketched in Fig. [10] (a) a realistic estimate of a part of the 
energy spectrum. The spectrum has mainly three con- 
stituents: (i) The thick (red) lines represent the eigenval- 
ues of the BHM of the first Bloch band (with U = U co "). 
The lower branch consists of three states with attracting 
particles (U < 0) in the same well, and the upper branch 
of three states with repelling particles (U > 0) in the 
same well. In the horizontal middle branch are all re- 
maining energies with particles in different wells. These 
states are only weakly influenced by the interaction, (ii) 
The bunch of thin lines going from top to bottom rep- 
resent "dimers" similar to the lowest BH states but in 
higher Bloch bands. Their energy consists of single- 
particle energies of Bloch states in different bands and the 
interaction energy [/ corr < 0. (iii) The thin horizontal 
lines at the top represent the energies of non-interacting 
states where one particle is in the first Bloch band and 
one in the second Bloch band of at least one direction. 
The light (blue) filling indicates that the energy regime 




FIG. 10: (Color online) Eigenenergies as a function of the scatter- 
ing length at ^ = 1.7. (a) Sketch of energy spectrum (see text 
for details), (b) Energy spectrum obtained from a full numerical 
calculation. The three thick (red) lines indicate energies which arc 
closest to the energies predicted by the BHM. 



of higher Bloch bands is reached. Note, for any depth of 
the lattice two repulsively interacting particles will enter 
this energy regime at the resonance as it is always half 
way to the energy of two particles in the second Bloch 
band. 

In Fig. [TU] (b) we have plotted in the region of reso- 
nance the CI eigenenergies of the lowest totally symmet- 
ric states in the trap. As they are of the same symmetry, 
avoided crossings occur. These appear close to where the 
estimated energy spectrum (Fig. [TO] (a)) shows (true) en- 
ergy crossings. The avoided crossings are caused by a 
more or less strong coupling between states in different 
Bloch bands, which is not incorporated in the estimated 
spectrum. Particularly, the energy of states which can be 
identified with BH states of repulsively interacting par- 
ticles is distorted due to avoided crossings. Thus, one 
cannot speak of pure repulsively interacting states in the 
first Bloch band any more. In contrast, the coupling 
between non-interacting BH states and the first exited 
"dimer" is very small and no energy distortion is visible 
at the resolution of Fig. [TO] (b). 

The complicated level structure indicates that the 
BHM, even if JJ C °" is used, has to be handled with care, 
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if repulsively interacting and also non-interacting states 
at a resonance are to be described. In its corrected form 
it is. however, useful as a first approximation to find level 
crossings. 

The rich structure of avoided crossings offers poten- 
tially a scheme to convert states of repulsively interact- 
ing, attractively interacting, and non-interacting nature 
into each other by performing sweeps of the scattering 
length. To estimate the dynamical behavior of the system 
the knowledge of the coupling strengths between states 
in different Bloch bands is crucial. Since these couplings 
are ignored within the BHM more sophisticated calcula- 
tions such as the one presented in this work have to be 
performed to gain a deeper insight. 



D. Extended Bose-Hubbard model 

Especially for shallow lattices both the interaction be- 
tween particles in nearest-neighbor (NN) wells and hop- 
ping to next-to-nearest neighbors (NtNN) can become 
important [13, HH, [HJ ■ One can easily account for this 
by including the according matrix elements from -ffwan of 
Eq. ([5]) into an extended Bose-Hubbard model (EBHM) 
with Hamiltonian 

Hebh =Hbh + J2{b i _ 1 k + b\b- 1 ) 

<tt> V ' (18) 

+ u-2 ]T (b^lkk + blblkhk) ■ 

<i,k> 

The NtNN hopping parameter J 2 = (w_i | £- + Vol \wi) 
decays much faster with the lattice depth than J. For 
Vq = hio it is one order smaller than J and for 
Vq = 1.5 hiu already two orders. Furthermore, the 
EBHM possesses in BH representation two new pa- 
rameters Ui = 47r/ ^ aBC J d 3 f w (r ) 3 wi(r ), and U 2 = 

^m" 5 " / d 3 rwo(r ) 2 wi{r ) 2 , describing NN interaction. 
These parameters are generally at least one order of mag- 
nitude smaller than U. Unfortunately, the BH repre- 
sentation of Ui and U2 cannot be improved like for U 
as there exists no analytical solution for two interacting 
particles in different wells. Therefore, only if the follow- 
ing conditions are met, the EBHM can be quantitatively 
better than the BHM: (i) The scattering length is small 
compared to the trap length (i.e. |a sc | <C Oho), since 
otherwise XJ\ and U2 are predicted inaccurately, (ii) The 
lattice is deep enough that errors of the BH parameters 
are smaller than U\, U2 or J 2 . 

Fulfilling the first condition, the case of weak interac- 
tion (a sc = — 176 ao) is discussed in the following. The 
second condition is not easy to meet in few-well sys- 
tems, since especially the BH representation of the onsite- 
cnergics reflects the optimal parameters only for deep lat- 
tices (see the inset in Fig. [5]), where U±, t/ 2 and J 2 are 



already two orders of magnitude smaller than U and J, 
respectively. 

However, many observables do not depend on the ab- 
solute value of the eigenenergies, but on their differences. 
Energy differences determine for example the character- 
istic frequencies of a state in a coherent superposition of 
different eigenstates. These frequencies can be measured 
experimentally with high accuracy such that deviations 
of the simple BHM can be observed [icl |. 

The influence of the errors of the onsite energies can be 
reduced by considering differences of eigenenergies as the 
latter do not depend on the absolute value of the on-site 
energies but only on the offset A = e\ — e§. Although A 
is also predicted inaccurately in the BH representation 
for Vo < 1-2 hw, it has often only a minor influence on 
the energy differences. This is especially the case for the 
energy difference E5 — E2 for which a diagonalization of 
the BH Hamiltonian without interaction yields 

E 5 - E 2 = VJ 2 + A 2 . (19) 

For the present potential with Vo < 1.2 Huj the relation 
J 2 ^> A 2 holds. Since the errors of the parameters J 
and U are small for lattice depths down to 0.7 fau, one 
can clearly identify effects of NN interaction and NtNN 
hopping as proper corrections of the BHM for this energy 
difference. 

V [E,] 

1 5 10 



T 




EBHM with V - V 

• • full calculation 



V [ho,] 

FIG. 11: (Color online) Energy difference E§ — E2 in different 
orders of approximation as a function of the lattice depth Vo for 
a sc fa — 180ao. The value of E5 — E2 is shown relative to the same 
energy difference predicted by the BHM in BH representation. The 
BH representation of to,t±\ was determined for koxg = 3ir/2 and 
the value of J and J2 for fcrj^B = 2-7T. 

In Fig. [TT] E$ — E2 is shown relative to its value for 
the BHM with parameters steming from the BH rep- 
resentation. The correct energy difference is compared 
to the one predicted by the BHM which is successively 
amended. First, only effects of NN interaction are incor- 
porated by the use of the EBHM without NtNN hopping 
(i.e. setting J 2 = 0). Then the full EBHM is considered 
which is finally corrected by replacing C/ BH by JJ cm . The 
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relative importance of these three corrections for differ- 
ent lattice depths becomes clearly visible. Due to the 
weakness of the interaction NtNN hopping and NN in- 
teraction have comparable influence, most dominant at 
Vq « 0.9 fojj. For deeper lattices the correction of the in- 
teraction parameter becomes increasingly important and 
dominates finally for Vb > 1.3 froj the correction of the 
BHM. 

By including all three corrections one is able to deter- 
mine the correct energy difference with about 1% accu- 
racy down to lattice depths of only Vq = 0.7 Hlj. Below 
this lattice depth the EBHM is unable to improve the 
BHM significantly. Additional effects steming from the 
confinement V CO nf finally dominate the properties of the 
energy spectrum and the description by a BHM or even 
an EBHM becomes inappropriate. 

Also for other differences of cigenencrgies the predic- 
tion of the EBHM improves the one of the BHM at least 
qualitatively (depending on the importance of A). Only 
for differences involving the first energy level this is not 
the case. Since for attractive interaction (as shown in 
Fig. II 1 [) particles in the ground state prcdominantely oc- 
cupy the middle well, the energy difference to higher lying 
states, where also the outhcr wells are occupied, is partic- 
ularly influenced by A. The error in the determination 
of A can therefore dominate over improvements of the 
EBHM. We have verified that the failure of the EBHM to 
improve the BHM results for energy differences involving 
the ground-state energy E\ does not occur, if repulsive 
interactions are considered. 



V. CONCLUSIONS 

We have investigated the validity of the standard and 
an extended Bose-Hubbard model in a triple- well optical 
lattice. We defined the optimal Bosc- Hubbard parameter 
set to be the one which reproduces best the eigenenergies 
of our full numerical calculations. This allowed to deter- 
mine when the Bose-Hubbard model is principally able 
to describe the system under consideration and whether 
corrections of the standard Bose-Hubbard parameters are 
needed. 

For shallow optical lattices the confinement to three 



wells had a notable impact on the hopping parameter 
and the onsite energies. This impact could not be fully 
reflected by the standard values of these parameters. Dis- 
agreements due to the confinement vanished, however, for 
lattice depths Vo larger than the ground-state energy of 
the harmonic approximation 1.5 hw. Therefore, in this 
regime our results are valid for both large optical lattices 
and few-well systems. For results concerning the inter- 
action parameter this is also true for shallower lattices, 
since this parameter is barely influenced by the confine- 
ment. 

The standard interaction parameter showed to be valid 
only when the scattering length is smaller than one per- 
cent of the trap length. However, for this parameter we 
found an analytic expression which approximates the op- 
timal interaction parameter very well for a large range of 
lattice depths and interaction strengths. Since its defini- 
tion is quite general, one can expect that it is extendable 
to hetero-nuclear systems, interacting particles in differ- 
ent Bloch bands, and lattices with different wavelengths 
in the three spacial directions. 

The error due to the neglect of the energy-dependence 
of the scattering length was quantified and found to be 
of less importance than effects caused by the incomplete- 
ness of the Wannier basis. For deep lattices this error is 
comparable to anharmonic effects. 

Furthermore, we found that the sextic expansion of 
one well can already cover the full onsite interaction in 
the lattice. We have presented and explained the avoided 
crossing of repulsively bound pairs with states in higher 
Bloch bands. Finally, an extended Bose-Hubbard model 
was investigated and it was shown that in a region of 
small lattice depths it can partially improve the standard 
Bose-Hubbard model significantly. 
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